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Oh 1 Abstract 

We study the production of photons through the non-equilibrium relax- 
ation of a disoriented chiral condensate. We propose that to search for non- 
equilibrium photons in the direct photon measurements of heavy-ion collisions 
^ ; can be a potential test of the formation of disoriented chiral condensates. 
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Recently, there have been many investigations into the formation of " Disoriented Chiral 
Condensates" (DCCs) following relativistic heavy ion collisions proposed by Bjorken et al. 
for a navel signature for the chiral phase transition The generally accepted picture 

of heavy energetic nuclei collisions is based on the Bjorken's scenario. In the collisions of 
these highly Lorentz contracted nuclei, they essentially pass through each other, leaving be- 
hind a hot plasma in the central rapidity region with large energy density corresponding to 
temperature above 200 MeV where the chiral symmetry is restored. This plasma then cools 
down via rapid hydrodynamic expansion through the chiral phase transition during which 
the long-wavelength fluctuations become unstable and grow due to the spinodal instabilities 
The growth of these unstable modes results in the formation of DCCs. The DCCs 
are the correlated regions of space-time where the chiral order parameter of QCD is chirally 
rotated from its usual orientation in isospin space. Subsequent relaxation of such DCCs to 
the true QCD vacuum is expected to radiate copious soft pions with a navel distribution in 
the ratio of neutral to charged pions, / (P(f) ~ 1/y/J), which could be a potential experi- 
mental signature of the chiral phase transition observable at RHIC and LHC. However, since 
these emitted pions will undergo the strong final interaction, the signals for the distribution 
P(f) may be severely masked and become indistinguishable from the background. It then 
becomes important to study other possible signatures of DCCs that would be less affected 
by the final state interaction. Electromagnetic probe such as photon and lepton with longer 
mean free path in the medium than the pions serves as a good candidate and can reveal 
more detailed non-equilibrium information on the DCCs with minimal distortion 

Minakawa and Muller || have recently suggested that the presence of strong electromag- 
netic fields in relativistic heavy ion collisions induces a quasi-instantaneous "kick" to the field 
configuration along the it direction such that it is plausible that the chiral order parameter 
in the DCC domains, if formed, will acquire a component in the direction of the neutral 
pion. In this Letter, we consider the production of photons through the non-equilibrium 
relaxation of a DCC within which the chiral order parameter initially has a non-vanishing 
expectation value along the 7r° direction and subsequently oscillates around the minimum of 
the effective potential. Our aim is to understand how the photons can be produced from this 
oscillating 7r° field via the dynamics of parametric amplification as well as spinodal instabili- 
ties. In Ref. ||, Boyanovsky et al. have extensively studied the photon production from the 
low energy coupling of the neutral pion to photon via the Ua(1) anomalous vertex. They 
have found that for large initial amplitudes of the 7r° field photon production is enhanced 
by parametric amplification. These processes are non-perturbative with a large contribution 
during the non-equilibrium stages of the evolution and result in a distinct distribution of 
the produced photons. Here we will take into account another dominant contribution that 
also involves the dynamics of 7r° due to the decay of the vector meson through the elec- 
tromagnetic vertex. Although the corresponding dimensionless effective coupling involving 
the vector meson is quite small perturbatively, as we will see later, in fact, for the large 
amplitude oscillations of the ir° mean field, the contribution to the photon production is of 
the same order of magnitude as the anomalous interaction. 

The relevant phenomenological Lagrangian density is given by 

C = Cfj + Cy + C^Oyy + Cy + Cvn^, (1) 

where 
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C a = ■ d$ - im 2 (t)$ • $ - A ($ • $) 2 + ha, (2) 

£ 7 + C^o^^i ~^ 22^2 ^ ^ ^apFnu, (3) 

£ y + £ V7r7 = - ^rnvVVp + ^-e a ^F aP (4) 

where i 7 )^ = c^A, — c^A^ is the electromagnetic field-strength tensor, and V^„ = d M V^ — c^V^ 
is the field-strength tensor of the vector meson with mass my. In addition, $ = (a, tt°, tt) is 
an 0(N+1) vector with n = (tt 1 , tt 2 , vr^ -1 ) representing the N — l pions. This Lagrangian 
without the vector meson piece is the model considered in Ref. ||. In this work, we will follow 
closely the formalism developed there. Likewise, we are going to ignore the hydrodynamical 
expansion and adopt the simple "quench" phase transition from an initial thermodynamic 
equilibrium state at a temperature (Tj) higher than the critical temperature (T c ) for the 
chiral phase transition cooled instantaneously to zero temperature. This "quench" scenario, 
which has been widely used in the study of non-equilibrium phenomena of DCCs |5|-f§, 
can capture the qualitative features of this non-equilibrium problem and allow a concrete 
analytical calculation. Thus, we take 0,03 



2 M m l 
m [t) = — 



2 



T{ > T c . (5) 



The parameters can be determined by the low-energy pion physics as follows: 

m a « 600 MeV, f n « 93 MeV, A « 4.5, T c « 200 MeV, 

h « (120 MeV) 3 , m v « 782 MeV, \ v ~ 0.36, (6) 

where V is identified as the u meson, and the coupling Ay is obtained from the uj — > 7r°7 
decay width |T0| . 

Since we are only interested in the photon production, we can integrate out the vec- 
tor meson to obtain the effective Lagrangian density that contains the relevant degrees of 
freedom given by 

2 \ 2 

C cS = C a + C 1 + - -i-^-e^e^ d x n%n°F^F af3 , (7) 

where the higher derivative terms are dropped out. At this point it must be noticed that here 
we have assumed the validity of the low-energy effective vertices. The effective vertices that 
account for the above mentioned processes may be modified in the strongly out of equilibrium 



situation [[LI]. In fact, one should obtain the non-equilibrium vertices by integrating out the 
quark fields and the vector meson in the context of the fully non-equilibrium formalism that 
we are currently studying in detail. 

Before proceeding further, let us discuss the qualitative features of the relative impor- 
tance of these two effective couplings to the photon production . For small amplitude 
oscillations of the 7r° mean field (e.g. smaller than the mass of the it ), from the naive 
perturbation argument with the parameters in Eq. (|6|), we expect that the effect of photon 
production from the coupling of (dir F) 2 is one order of magnitude smaller than that of the 
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anomalous interaction. However, for large amplitude oscillations, the non-linearity results in 
the (dn°F) 2 term being of the same order of magnitude as the anomalous interaction, and 
moreover this non-linear effect will make the photon production mechanism due to these 
two couplings more effective as we can see later. 

Following the non-equilibrium quantum field theory that requires a path integral repre- 
sentation along the complex contour in time |12|], the non-equilibrium Lagrangian density is 
given by 



£ neq = £ eff [$ + ,A+] - C eS [$-,A 



(8) 



where +(— ) denotes the forward (backward) time branches. The non-equilibrium equations 
of motion are obtained via the tadpole method. As mentioned before, the situation of interest 
to us is a DCC in which both the a and 7r° fields acquire the vacuum expectation values. 
We then shift a and ir° by their expectation values described by the initial non-equilibrium 
states specified later, 



a(x, t) = (f)(t) + x{x, t), (j>(t) = (a(x, £)), 

with the tadpole conditions, 

( x (x,t)) = 0, &(x,t)) = 0, (vf(f,t))=0. 



(9) 



(10) 



This tadpole conditions will be imposed to all orders in the corresponding expansion. In order 
to derive the non-equilibrium evolution equations that incorporate quantum fluctuation 
effects from the strong a— ir interactions, we will use the large- N limit to provide a consistent, 
non-perturbative framework to study this dynamics. To leading order in the 1/N expansion, 
following the Hartree factorizations (Eqs. (2.7)-(2.9) of Ref. [§]) implemented for both ± 
components, the Lagrangian then becomes 



C 



eil' 



c 



cff 



+ x X + 4> ,7? ,A: 



•^M 2 (t) X + 2 - \Ml{t)^ 



-Ml(t)n+ 2 - -F+F+^ + 



32tt 2 A 



2 

e 2 X 2 



v 



iC(t)) 2 e^ 5 e^ s F+F+ + 
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V 



e 2 X 



4mym 2 



e 2 X 2 



v 



4mym 2 



at)^ os e a ^ d a F+F ap 



e 2 X 2 



m v m< 



_#»»ffh d^+d^+F+F^ - {+ -> -} 



11^ 



where 



Ui(t) = 4>{t) + [m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + 4A£(t) <f)(t) - h 
U 2 (t) = C(t) + \m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + 4AE(t)l 0), 
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M 2 x (t) 



m 2 (t) + 12A0 2 (t) + 4AC 2 (t) + 4A£(t) 
AfJ(i) = m 2 (t) + 4A0 2 (t) + 12AC 2 (t) + 4A£(t), 
MH(t) = m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + 4A£(t), 
S(t) = (vr 2 )(t)-( 7 r)(0). 



(12) 



The expectation value described by the initial non-equilibrium states will be determined self- 
consistently. Here, we have performed a subtraction of (vr 2 )(t) at t = absorbing (7? 2 )(0) 
into the finite renormalization of the mass term. 

Since we consider the direct photon production driven by the time dependent oscillating 
field in which the photon does not appear in the intermediate states. It proves to be 
convenient to choose the Coulomb gauge that contains physical degrees of freedom without 
any other redundant fields f|. With the above Hartree-factorized Lagrangian in the Coulomb 
gauge, following the tadpole conditions, we can obtain the full one-loop equations of motion 
while we treat the weak electromagnetic coupling perturbatively: 



4>{t) + [m 2 {t) + 4A0 2 (t) + 4AC 2 (t) + 4A£(t)] <f>{t) 

C(t) + [m 2 (t) + 4A0 2 (t) + 4AC 2 (t) + 4A£(t)] Q{t) - 

e 2 \ 2 v d 



h = 0, 

„2 



c a(3fiv 



m 2 my dt 



C(t)(^V 2 A' T )(t) 



0. 



(13) 



Now we decompose the fields n and A? into their Fourier mode functions UAt) and 



VyAt) respectively, 



7f(x, t) 

A T {x,t) 



d 3 k 



A=l,2 



2(2vr)3^ 
d 3 k e 



%^(t)e^ + h.c 



(14) 



Ak 



where and b x j: are destruction operators, and e A £ are circular polarization unit vectors. 
The frequencies w z and uJ A t can be determined from the initial states and will be specified 
below. Then the mode equations can be read off from the quadratic part of the Lagrangian 
in the form 



dt 2 



+ k 2 + m 2 {t) + 4A0 2 (t) + 4AC 2 (t) + 4A£(t) 



e 2 A 2 



1 - 



dt 2 



V 2k {t) + 



1 - 



e 2 A 2 



m 2 my 



mtmt 



C 2 (t) 



C 2 (t) 



k 2 V lk {t) 



2n 2 U 



U k (t) 

C(t)v lk (t) 



0. 



k 2 V 2k {t) + k-—((t)V 2k (t) = 0. 



(15) 



With Ay = 0, this reproduces the mode equations derived in Ref. ||. To solve these mode 
functions, we must specify initial conditions. At the time of "quench", we assume that 
the quantum fluctuations for the 7r fields which undergo the strong interactions are in the 



5 



local thermodynamic equilibrium at the initial temperature T{ > T c with the chiral order 
parameter displaced initially away from the equilibrium position and with a non-vanishing 
initial amplitude along the n° direction, i.e., £(0) 7^ and 0(0) = 0. However, since the 
photons interact electromagnetically, their mean free paths are longer than the estimated 
size of the presumed quark-gluon plasma fireball so that the produced photons will escape 
from the plasma freely Therefore, one can argue that the photonic medium effects play 
no role in the dynamics of photon production. Based on the above arguments, the initial 
conditions for the mode functions are given by 



£4(0) = 

V Xk (0) 



17* (0) = 
, Vxk(0) 



-iu„ k , u4 = k 2 + m 2 (t < 0) + 4A[0 2 (O) + C 2 (0)] 
—iu Ak , u Ak = k, 



(16) 



with the expectation values with respect to the initial states given by 



E(t) = (N - 1 

fA 



A d 3 k 



\u k (t)\ 2 -i 



coth 



2(27r) 3 ^ fc 

d 3 k r\ r 

" \V2k(t)\ 2 - \V lk (t)\ 2 



(^V 2 A J T )(£) 



A d 3 k 
2(2tt)W 



dt 



\V lk (t)\ 2 + \V 2k (t)\'< 



(17) 



where we set the cutoff scale A ~ my and N = 3. The above specified initial conditions are 
physically plausible and simple enough for us to investigate a quantitative description of the 
dynamics. The expectation value of the number operator for the asymptotic photons with 
momentum k is given by [H 



(N*(t)> 



1 

2k 



A T (k, t) ■ A T (-k, t) +k 2 A T (k, t) ■ A T (-k, t) 



- 1 



^E[l^(t)| 2 + A; 2 |F Afc (t)p 



- 1. 



(18) 



This gives the spectral number density of the photons produced at time t, dN(t)/d 3 k. 

We now perform the numerical analysis. We choose to represent the quench from an 
initial temperature set to be Tj = 220 MeV = 1.1 fm -1 to zero temperature ||. The 
evolution is tracked up to a time of about 10 fm after which the hydrodynamical expansion 
becomes important. Fig. |1] shows the temporal evolution of the ((t) by choosing the initial 
conditions C(0) = 0.5 fm -1 and 1 fm" 1 , and C(0) = (f)(0) = 0(0) = 0. The ((t) evolves 
with damping due to the backreaction effects from the quantum fluctuations. In Fig. |2], we 
present the time-averaged invariant photon production rate, kdR/d 3 k, where 



dR=- 
T 



dN(t) 
dt 



dt, 



(19) 



over a period from the initial time to time T = 10 fm. 

In the case of £(0) = 1 fm -1 , the produced non-thermal photons have spectrum peaks 
around two photon momenta, k = 1.35 fm -1 and k = 2.7 fm -1 , which exhibits the features 
of the unstable bands and the growth of the fluctuating modes. The growth of the modes 
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in the unstable bands translates into the profuse particle production. Note that the peaks 
are located at k = and k = uo^ where uo^ is the oscillating frequency of the ((t) field 
with £(0) = 1 fm" 1 in Fig. [D. The spectrum peaks clearly result from the oscillations of the 
( field that serves as the time dependent frequency term in the mode equations of At (fL5"l) . 
Thus, the photon production mechanism is that of parametric amplification. Comparing 
with the results in Ref. ||, where the authors consider the photon production only via the 
C/a(1) anomalous interaction, and the dotted curve in Fig. |2| which denotes the time-averaged 
invariant photon production rate for £(0) = 1 fm -1 with the vector meson channel turned 
off (Ay = 0), we can easily recognize that the 1.35 fm -1 peak is resulted from the coupling 
w°FF while the 2.7 fm -1 peak is from the interaction (dn°F) 2 . For a smaller initial field 
amplitude £(0) = 0.5 fm -1 , the oscillating frequency decreases while the peaks shift to the 
lower-momentum region with peak values almost two orders of magnitude lower than those 
of C(0) = 1 fm" . This explicitly shows the non-linearity of the amplification process. 

We now provide the analytical analysis to understand qualitatively the above numerical 
results and especially why the (dir°F) 2 coupling is dominant at k — 2.7 fm -1 in photon 
production although it is small perturbatively. From Fig. [I] for £(0) = 1 fm -1 , it shows that 
the solution of £(t) is a quasiperiodic function with a decreasing amplitude during the first 
few oscillations. To obtain the analytical estimates for the locations of the unstable bands 
in the produced photon spectrum as well as their growth rates, let us approximate the ((t) 
as ((t) ~ C 5111(0;^^) . The ( is the average amplitude over a period from the initial time up 
to time of 10 fm, and is about ( ~ 0.8 fm -1 , and the oscillation frequency, ~ 2.7 fm" 1 , 
measured directly from Fig. |l|. Then, the photon mode equation in Eq. fll5[) becomes 



2 



2m 2 m 2 , 



cos 



(2u; c t) 



'■v 



k 2 V lk {t) - kp^f cos{u c t)V lk {t) = 0. (20) 



When the vector meson channel is turned off (Ay = 0), we change the variable to z = u^t/2. 
Then, Eq. (^) becomes 

d 2 Ak 2 Ake 2 C 

i? v ^^~2^-r a[2z)Vit=0 - (21) 

This is the standard Mathieu equation |14]]. The widest and most important instability 



is the first parametric resonance that occurs at k = lo^/2 with a narrow bandwidth 5 ~ 
e 2 (/ (2ir 2 f n ). The instability leads to the exponential growth of photon modes with a growth 
factor / = e 2)iz , where the growth index fi ~ 5/2. This growth explains the peak at 
k = 1.35 fm -1 in Fig. |2|. When the Ua(1) anomalous vertex is turned off (f w — > 00), we 
change the variable to z' = uj(t. Then, Eq. ( p0|) becomes 

if 2 k 2 k 2 p 2 \ 2 ? 2 

V lk + Vi* - 22 cob(20 V« = 0. (22) 



dz' 2 uj 2 2m 2 my 



Now, the parametric resonance occurs at k = uo^ with a growth factor /' = e 2 ^' 2 ', where 
// ~ e 2 \y( 2 uj 2 / (8mjmy). This growth explains the peak at k = 2.7 fm -1 in Fig. 0. Taking 

( = 0.8 fm -1 and t ~ 10 fm, the ratio of their growth rates is given by /'// — 0.5. This 
means that the height of the 2.7 fm -1 peak is about one half of the 1.35 fm" 1 peak as we 
can see in Fig. [2|. 
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We then compare our results with the thermal photon emitted from a quark-gluon plasma 
and a hadron gas. In Fig. |2|, the invariant photon production rate for the quark-gluon plasma 
is drawn using the parameters given in Ref . . For the hadron gas, we have used the rates 



for the most important scattering and decay processes flql - It is shown that the 1.35 fm~ 
peak and the 2.7 fm~ peak is about an order of magnitude larger than the thermal photons. 
Therefore, we can come to the conclusion that these non-thermal photons can be regarded 
clS db distinct signature of non-equilibrium DCCs. 

In conclusion, we have studied the production of photons through the non-equilibrium 
relaxation of a disoriented chiral condensate within which the chiral order parameter initially 
has a non-vanishing expectation value along the 7r° direction . Under the "quench" approxi- 
mation, the invariant production rate for non-equilibrium photons driven by the oscillation of 
the 7T° field due to parametric amplification is given, which exceeds that for thermal photons 
from a thermal quark-gluon plasma or hadron gas for photon energies around 0.2 — 0.7 GeV. 
These relatively high-energy non-thermal photons can be a potential test of the formation 
of disoriented chiral condensates in relativistic heavy-ion-collision experiments. 
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FIGURE CAPTIONS 



Fig. |T] Time evolution of the mean field £(£) with initial amplitudes 0.5 fm 1 and 1 fm 1 . 

Fig. |2] Lower and upper solid curves are the time-averaged invariant photon production rates 
for C(0) = 0.5 fm -1 and 1 fm -1 respectively. The latter shows profuse photon production. 
The dot-dashed curve is the invariant photon production rate for the thermal hadron gas 
at T = 220 MeV. The rate for T = 220 MeV quark-gluon plasma is denoted by the dashed 
curve. Also shown is the dotted curve which denotes the rate for £(0) = 1 fm~ with the 
vector meson channel turned off, i.e., Xy = 0. 
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